Equilibria for a gyrostat in newtonian interaction with 

two rigid bodies 



J. A. Vera 
Universidad Politecnica de Cartagena 
Departamento de Matematica Aplicada y Estadfstica 
Paseo Alfonso XIII, 52, 30203 Cartagena (Murcia), Spain 
e-mail: juanantonio.vera@upct.es 

Abstract 

In this paper the non-canonical Hamiltonian dynamics of a gyrostat in the three 
body problem will be examined. By means of geometric-mechanics methods some 
relative equilibria of the dynamics of a gyrostat in Newtonian interaction with two 
rigid bodies will be studied. Taking advantage of the results obtained in previous 
papers, working on the reduced problem, the bifurcations of these relative equilibria 
will be studied. The instability of Eulerian relative equilibria if the gyrostat is close 
to a sphere is proven. Necessary and sufficient conditions will be provided for 
lineal stability of Lagrangian relative equilibria if the gyrostat is close to a sphere. 
The analysis is done in vectorial form avoiding the use of canonical variables and 
the tedious expressions associated with them. In this way, the classic results on 
equilibria of the three-body problem, many of them obtained by other authors who 
had used of more classic techniques, are generalized. 

AMS Subject Classification: 34J15, 34J20, 53D17, 70F07, 70K42, 70H14 

1 Introduction 

In the last few years a lot of papers about the problem of roto-translational motion of 
celestial bodies have appeared. These show a new interest in the study of configurations 
of relative equilibria by differential geometry methods instead of more classical ones. 

In the present paper, references will be made to Wang et al. [H] , concernig the problem 
of a rigid body in a central Newtonian field, and Maciejewski concerning the problem 
of two rigid bodies in mutual Newtonian attraction. Let us remember that a gyrostat is 
a mechanical system S composed of a rigid body <S" and other bodies S" (deformable or 
rigid), connected to it in such a way that their relative motion with respect to its rigid 
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part do not change the distribution of mass of the total system S (Leimanis, j2]). These 
papers have been generalized by Mondejar and Vigueras [4J to the case of two gyrostats 
in mutual Newtonian attraction. 

With regards to the problem of three rigid bodies, references will have to be made to 
Vidiakin 8j and Duboshin pQ who proved the existence of Euler and Lagrange configura- 
tions of equilibria when the bodies possess symmetries. Later Zhuravlev and Petrutskii 
[TU] reviewed the results until 1990. 

Vera |Hj and Vera and Vigueras (|E],[Zj) have studied the non-canonical Hamiltonian 
dynamics of n + 1 bodies in Newtonian attraction, in which n of them are rigid bodies 
with a spherical distribution of masses or material points and the other one is a triaxial 
gyrostat. Working in the reduced problem, global considerations about the conditions 
for relative equilibria were made. Finally, in an approximated model of the dynamics, a 
study of some relative equilibria of a gyrostat in Newtonian interaction with three rigid 
bodies is carried out. 

In this paper, the existence and stability of some periodic solutions of the dynamics of 
a gyrostat in Newtonian interaction with two rigid bodies will be investigated. In a first 
approach to the problem, the attitude dynamics of both rigid bodies are assumed to be 
the same as the dynamics of a rigid body in torque free motion. The two rigid bodies have 
revolution symmetry about the third axis of inertia. On the other hand, the Newtonian 
interaction of the gyrostat between both rigid bodies is simplified supposing the gyrostat 
a material particle. The Newtonian interaction between the two rigid bodies are the same 
of two material particles. 

With these hypotheses, according to Vera and Vigueras [Zj, working in the double 
reduced space of configuration of the problem, the equations of motion and those which 
determine the relative equilibria will be derived. 

Two families of relative equilibria, Eulerian and Lagrangian will be investigated. The 
Eulerian relative equilibria are completely determinated by a polynomial equation of de- 
gree nine. The bifurcations of these equilibria are made. The bifurcations of the La- 
grangian relative equilibria are completely investigated and expressions for the Lagrangian 
relative equilibria when both solids are close to spheres are provided. 

We should notice that the studied system, has potential interest both in astrodynamics 
(dealing with spacecrafts) as well as in the understanding of the evolution of planetary 
systems recently found (and more to appear), where some of the planets may be modeled 
like a rigid body rather than a rigid body. In fact, the equilibria reported might well be 
compared with the ones taken for the 'parking areas' of the space missions (GENESIS, 
SOHO, DARWIN, etc) around the Eulerian points of the Sun-Earth and the Earth-Moon 
systems. Some interesting numeric results have been calculated in the Appendix B. 

The methods used in this work can be applied to similar problems. The study of the 
nonlinear stability of the relative equilibria obtained here is the logical continuation of 
this work. 
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2 Poisson dynamics 

Let Sq be a gyrostat of mass m , and Si, S2 two symmetrical rigid bodies of masses mi 
and m 2 respectively; 3 = {O, ui, u 2 , u 3 } an inertial reference frame; 3 = {Co, bi, b 2 , b 3 } 
a body frame fixed at the center of mass C of S , (see Vera and Vigueras [7] for details). 

The following notation is used 

. , mima m Q M 2 

M 2 =m 1 + m2, Mi = mi + m 2 + m Q , #1 = ^7 — , 92 



M 2 ' ^ Mi 

Where for u, v e M 3 , u ■ v is the dot product; | u | is the Euclidean norm of the vector 
u; u x v is the cross product; Ir3 is the identity matrix; and is the zero matrix of order 
three. Consider I = diag(A , Aq, Co) the diagonal tensor of inertia of the gyrostat, and 
Ij = diag(Ai, A*, CA the diagonal tensors of inertia for the rigid bodies Si, i = 1, 2. The 
generic expression z =(n~i, n 2 , n , A, p\, /x, p M ) e R 21 is an vector of the twice reduced 
problem obtained by appling the symmetries of the system. n = lo^o + is the total 
rotational angular momentum vector of the gyrostat in the body frame, which is attached 
to its rigid part and whose axes have the direction of the principal axes of inertia of Sq; 
the vector \ r = (0, 0, 1) is the constant gyrostatic momentum and 11, = Ijfij (i = 1, 2) are 
the total rotational angular momentum vectors for the two rigid bodies. The elements 
A, [i, p\ and p M are respectively the barycentric coordinates and the linear momenta 
expressed in the body frame 3- 

Following the results of Vera and Vigueras [Ij , according to the hypotheses formulated 
in the introduction of this paper, a good approximation to the potential of the system is 
expressed by the following expression 



where 



V = Vi+V 2 



. Gmim 2 Gmimo Gm 2 mo 

Vl = - , x 1 1" 1 mo x 1 + 



... . Gm «i Gm a 2 3Gm /i 3Gmof 2 

V2 = — - 1 ^„ > 10 + 



M-^A|3 |/x+^A|3 \n-vg\\?> 1^+1^ I 5 



and 



ai = 2Ai + C u a 2 = 2A 2 + C- 
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The Hamiltonian function of the system adopts the form 



H(z) = + + in i % - i r • i x n + Jnjr'n! + Vi 2 ~ x n 2 + v(a, /*) 

2(yfi Z# 2 Z 11 

Let (M, { , },H) with M = IR 21 be the Poisson manifold, where { , } 
is the Poisson brackets defined by means of the Poisson tensor 



B(z) 
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In B(z), v is considered to be the image of the vector v G 
phism between the Lie Algebras R 3 and 30(3), i.e. 

-v 3 v 2 
v = I f 3 —vi 
-v 2 v 1 



by the standard isomor- 



The equations of the motion are given by the below expression 

{z,H(z)} = B(z)V z H(z) 



dz 
~dl 



with V Z V being the gradient of V with respect to an arbitrary vector z. 

Calculating {z,7Y(z)}, the below group of vectorial equations of the motion can be 
written as 

° 1 1„ x Q • A x V A V • //. x V /( V 



dt 

^A pa » „ dp\ ^ _ , , 

— = — + \xn , -f = PA xn„-v A v 

dt gi dt 



dfi _ p M 



dt 

dUi 



+ fl x fl 



dpfj. 



(1) 



92 

-- n x x o x , 



°' dt 

dU 2 



= p^o - v M v 
n 2 x n 2 



dt ' ' dt 

Important elements of B(z) are the associated Casimir functions. The vector 



L = II + Axp A + /jxp )1 
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is a part of the total angular momentum L given by 

L = n 2 + IIx + L 
Then the below result can be concluded. 

it i2 in ■ I 2 

Proposition 1 If (p^ (i = 0, 1, 2) are real smooth functions, then <£>o( 2~)j fi\2~)i 
(i = 1, 2) are Casimir functions of the Poisson tensor B(z). Furthermore, Ker~B(z) =< 
V z <£> , 'V z (p l , V z <£> 2 >. We also have %■ = 0, which means the total angular momentum 
vector remains constant. If IIq = (7rJ, 7Tq, tTq), then 7Tq is an integral of the motion. 



3 Relative Equilibria 

If z e = (II|, Ilf, IIq, A 6 , p A , /x e , p£) is a generic relative equilibrium, the below 
vectorial equations are verified 

ITq x Sl e + A e X (V A V) e + /x x (V M V) e = (2) 



^ + A e xOS = 0, pi x n§ = (V A V) e 

^ + M e x = o, P « x n« = (v M v) e 

#2 

iTi x n\ = o, x n e 2 = o (3) 

where (V A V) e and (V M V) e are the values of V A V and V M V in z e . 

According to the relationships provided by Vera and Vigueras [Zj, the following results 
are obtained. 

Lemma 2 Whenever z e = (II|, ITj, IIq, A 6 , p A , n e , p^) zs a relative equilibrium, the 
below relationships are verified 

I ng | 2 | a £ | 2 -(A e • n^) 2 = -(A e • (v A v) e ) 

|^| 2 |/x e | 2 -( M e -^) 2 = i(/i e -(V M V) e ) 

#2 

The previous two identities will be used to obtain necessary conditions for the existence 
of relative equilibria. 

Certain relative equilibria will be studied assuming that vectors fig, A e and /z e satisfy 
special geometric properties. 

Definition 3 z e is said to be an Eulerian relative equilibrium when, A 6 and fi e are 

proportional and Q e is perpendicular to the straight line that is generate. 
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Definition 4 z e is said to be a Lagrangian relative equilibrium if X e and fi e are not 

proportional and Q e is perpendicular to the plane that both of them generate. 

From the above definitions, the following property is deduced. 

Proposition 5 In an Eulerian or Lagrangian relative equilibrium, momenta are not ex- 
ercised over the gyrostat. 

Next, necessary and sufficient conditions for the existence of Eulerian and Lagrangian 
relative equilibria will we obtained. 

4 Eulerian relative equilibria 

According to the relative position of the gyrostat So with respect to Si and S 2 , there 
are three possible equilibrium configurations: a) S 2 SiS , b) S 2 S Si and c) S S 2 Si. 

4.1 Necessary conditions of existence 

The following lemma is a direct consequence of the geometry of the problem. 

Lemma 6 If z e = (II!, Ilf, LTq, A e , p^, /x e , p^) is a relative equilibrium of Euler type, 
then for the configuration SqS 2 S\ 

fi ——j-X = A e + /j, +—j-\ 

1 M 2 11 1 1 M 2 1 



In a similar way, for the configuration S 2 S S 



1 1 1 M 2 1 1 M 2 
Finally, for the configuration S 2 SiSo 

1 M 2 1 1 M 2 
If z e is an Eulerian relative equilibrium, them 

9l \fl° e | 2 |A e | 2 =A e -(V A V) e 

»0 1 2 1 ..e |2 



92 \nl |V e |^ M e .(V M V) f 



with 



u e_ m 2 x e _ ye e , m lu /i , fl]) .e e _ (i 1 + P) m 2 + » 
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where p G (— oo, — 1) in the configuration a), p G (—1,0) in the configuration b) and 
p G (0, +oo) in the configuration c). Moreover the following expressions are possible to 
be obtained 

(V A V) e = h l {p)X e , (V M V) e = h 2 (p)X e 



with 



h(p) 



Gm\vri2 Gm§m\sgn{\ + p) 



I * e I 3 
Gmom2sgn(p) 

Mo I A e I 3 



M 2 , 
mi 

~7 + 7 



\ e i 3 

Pi 



m 2 



+ 



1 + p) 2 (1 + p) 4 | A 



e 12 



A e I 2 



and 



Hp) 



Gm sgn(l + p) 



Gm sgn(p) 
I A e I 3 



m 2 



02 



mi 

— + 

P 2 



(1 + P) 2 

gi 
A e | 2 p 4 



i e 12 



+ 



(4) 
(5) 

(6) 
(7) 



where ^ = 3(Ci - Ai)/2, /3 2 = 3(C 2 - A 2 )/2 and san(x) 
Now, from the identities 



A e • (V A V) e 
» e ■ (V M V) £ 

these equations are deduced 



n e i 2 - 



r> e i 2 - 



A e ^(p) 

((1 + p)m 2 + pmi) 



,e |2 



Mo 



1 if x > 
-1 if x < 



mi + m 2 



771x1712 

77l + 77l\ + m 2 



-^2(P) 



m ((1 + p)m 2 + pmi) 

Then for an Eulerian relative equilibrium, p must be a real root of the below equation 

m (mi + m 2 ) ((1 + p)mi + pm 2 ) /ii(p) = mim 2 (m + mi + m 2 )h 2 (p) (8) 

The following proposition summarizes all these results. 

Proposition 7 If z e = Ilf, IIq, A e , p A , /x e , p^) zs an Eulerian relative equilibrium, 
the equation (0) /ias, a£ /east, a rea/ root in which the functions h\{p) and h 2 (p) are given 
by the expressions Ulfy). The modulus of the angular velocity of the gyrostat is 

rrii + m 2 



n e i 2 - 



'MP) 



Remark 8 When | A e | has a fixed value, if an Eulerian relative equilibrium exists, the 
equation (Q5J) has real solutions. The number of real roots of the equation (Q|) will depend, 
obviously, on the parameters which exist in the system. 



7 



J. A. Vera 



Equilibria for a gyrostat in newtonian interaction with two rigid bodies 



4.2 Sufficient conditions of existence 

The below proposition indicates how to calculate solutions of eq. (j2J). 

Proposition 9 When | A e | has a fixed value, let p be a solution of the equation (0), 
where the functions hi(p) and h 2 (p) are given by the expressions yi6|). Then z e = (Ilf, 
nf, ng, A e , p|, /x e , p«) #wera 6y 

A e = (A e , 0, 0), M e = (/A 0, 0), fig = (0, 0, wg) 

P ^ = (o,^A e ,o), P ^ = (o,^/i e ,o), ng = (o,o,c ^ + o 



,e 



((1 + p)m x + pm 2 ) 2 _ (mi + m 2 )h 1 (p) 



v = u A ' ^ 

M 2 mim 2 
is an Eulerian relative equilibrium. The total angular momentum of the system is ex- 
pressed by 

L = (0, 0, C 2 u e 2 + C x u\ + C ul + I + gi uj e (X e ) 2 + g 2 u e (p e ) 2 ) 
with I being the gyrostatic momentum. The vectors Yl\ verify the vectorial equations 

n;xn; = o, n^xfi; = o 

4.3 Eulerian relative equilibria when S2 and S\ are spherical 
rigid bodies 

Consider the existence and number of solutions for Eulerian relative equilibria when 
62 and Si are spherical rigid bodies. In this case C x = A\, C 2 = A 2 and the equation (jHJ) 
is equivalent to the below polynomial equation 

(mi + m 2 )p 5 + (3m 2 + 2nii)p A + (3m 2 + m x + m (sgn(l + p) (9) 
—sgn(p)))p 3 + (777,2 — m 2 sgn(l + p) — misgn(p) — 3m sgn(p))p 2 

— (3sgn(p)niQ + 2sgn(p)nii)p — (m sgn(p) + m\sgn(p)) = (10) 

This equation has a unique real solution in the intervals (— 00, — 1), (—1,0) and 
(0, +00). Therefore only one Eulerian relative equilibrium exists. 
On the other hand 

2= Gjnii + m 2 ) / m ( sgnjl + p) _ sgn(p) \ 

1 1 I A e | 3 V (mi + m 2 ) V (1 + p) 2 P 2 J 

where p is the only solution of the equation 0. 

Proposition 5 gathers the results about Eulerian relative equilibria when S 2 and Si 
are spherical rigid bodies in the configurations a), b) and c). 
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Proposition 10 1. If p is the unique positive root of the equation 

(mi + m 2 )p 5 + (3m 2 + 2m 1 )p 4 + (3m 2 + m 1 )p 3 + 
— (3m + mi)p 2 — (3m + 2m\)p — (m + mi) = 

with 

|2 _ G(m l +m 2 ) , m f \ 



"o 



A e | 3 V K + m 2 ) V(l + P) 2 P 2 
i/ien z e = (n|, n^, n§, A e , p A , /x e , p*), </u>en fry 

A e = (A e ,o,o), M e = (p e ,o,o), ns = (o,o,wS) 
P e A = (o, ^A e , o), P ^ = (o, o), ng = (o, o, c u e + 1) 

is the unique solution of relative equilibrium of Euler type in the configuration 
5*2 5*1 5*0 ■ 

2. If p G (—1,0) is the unique root of the equation 

(m 1 + m 2 )p 5 + (3m 2 + 2m 1 )p 4 + (3m 2 + m 1 )p 3 + 
+ (3m + 2m 2 + mi)p 2 + (3m + 2rai)p + (m + mi) = 

with 

e 2 _ G(m! + m 2 ) ^ m ( 1 1 



1 01 ~ |A e | 3 V K + m 2 ) V(l + p) 2 p 2 
then z e = (II|, Ilf , ng, A e , p A , /x e , p£) #wen fry 

A e = (A e ,o,o), M e = (p e ,o,o), n e = (o,o,u e ) 

p e x = (0, ^A e , 0), p^ = (0, g 2 u e p e , 0), W = (0, 0, C u e + /) 

£/ie unique solution of relative equilibrium of Euler type in the configuration 
S2S0S1. 

3. If p e (—00, —1) £/ie unique root of the equation 

(mi + m 2 )p 5 + (3m x + 2m 2 )p 4 + (2m + 3m 2 + m 1 )p 3 + 
+ (3m + mi)p 2 + (3m + 2m 1 )p + (m + m x ) = 

where 

2 G(mi + m 2 ) ^ ; m 



1 A e I 3 V K + m 2 ) VP 2 (1 + p) 2 , 
toen z e = (II|, Ilf, ng, A e , p A , /x e , pp #wen fry 

A e = (A e ,0,0), M e = (p e ,0,0), ^ = (0,0,^) 
p e A = (0, 9l u e X e , 0), p^ = (0, g 2 u e Q p e , 0), Ilg = (0, 0, C u e + Z) 

£/ie unique solution of relative equilibrium of Euler type in the configuration 
S0S2S1. 
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4.4 Eulerian relative equilibria when S2 and Si are not spherical 
rigid bodies 

In the present case, after carrying out the appropriate calculations, the equation (JHJ) 
is reduced to the study of the positive real roots of the nine degree equation 

(3 2 q{p) - m x m 2 a 2 p 2 {p + l) 2 p(p) = (11) 

where 

p(p) = (mi + m 2 )p 5 + (3m 2 + 2m 1 )p A + (3m 2 + m 1 + m (sgn(l + p) (12) 
—sgn(p)))p 3 + (777,2 — 7772.8(777(1 + p) — misgn(p) — 3mosgn(p))p 2 

— (3sgn(p)rriQ + 2sgn(p)mi)p — (m sgn(p) + m 1 sgn(p)) (13) 

and 

q(p) = mo(ksgn(p)m 2 — sgn(l + p)mi)p 5 + m 2 (sgn(l + p)777! + 5sgn(p)mok + kmisgn(p))p 4 
+2ksgn(p)(2mi + 5mo)p 3 + 2ksgn(p)(3mi + hmo)p 2 + ksgn(p)rri2( / kmi + 5m )p 

+ksgn(p) 7772(7770 + 777^ 

where 

(3, = 3(d - AO/2, (3 2 = 3(C 2 - A 2 )/2 

with /3j = fc/3 2 , a =| A e I and fcel. 

In order to study the number of real roots of the polynomial ([TTj) the rational function 
will be studied 

m 1 m 2 a 2 p 2 (p + l) 2 p{p) 

p2 = {p) = W) 

In practical applications 777 is very small, then up to the first order in 777 
o p / x aV(p + l) 2 Pi(p) , v 

where 

pi(p) = p 5 + (2 + p)p 4 + (1 + 2p)p 3 + (p - p^r7(l + p) 
-(1 - p)sgn(p))p 2 - 2sgn(p)(l - p)p - (1 - p)sgn(p) 

and 

Qi(p) = { s 9 n {l + p) + s 9 n (p)k)p A + 4:sgn(p)kp 3 + 6sgn(p)kp 2 + Asgn(p)kp + sgn(p)k 
where p = . 

777 2 + 777! 
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The polynomial qi has no roots in (0, +00) and (— 1, 0) if k > and k < 0, respectively. 
On the other hand, qi has only one root, p l , in (0, +00) and (—1, 0) if k < and k > 0, 
respectively. p will be denoted as the only root of p±. The implicit curve Res(k,fi) = 0, 
where Res is the resultant of the polynomials p± and qi, is used to study the graph of R±. 
When /z has a fixed value, the only /c which verifies Res(ko, fi ) = 0, exists according to 
the Implicit Function Theorem. Consider the only root, p\ = p (/i ), of qi for k and fi . 
The expressions 

(Pmax> fl(*0 = #l(/W))> (Pmin, &(*0 = ^l(Pmin)) 

are the local maximum and minimum of the function Ri. 
In the configuration S 2 SiS , for any value of /i fixed 

lim f 2 (A;) = 0, lim . £ 2 (h) = -00 



if k > 0. For k < 



lim £i(A;) = 0, lim £ 2 (&) = +°° 

fc-++oo A;— >0+ 



if Pi > Pi- If Pi < Pi, R\ is strictly increasing. 

For the configuration S 2 S Si if < p x and A; > 0, the function Ri has just a minimum, 
which verifies 



If Pi < Pi, then 



If k < 0, then i?i verifies that 



lim £ 2 (/c) = £ 



lim f 2 (]fc) = 

fc— >+oo 



lim £ 2 (k) = £ , lim ^(fc) = +00 
lim = 0, lim ^i(fc) = 

The results for the configuration S S 2 Si will be deduced from the configuration S 2 SiS . 
According to these statements, the following proposition can be stated. 



Proposition 11 In the configuration S 2 SiSq the following results are verified. 
a) For k > then: 

1. If (3 2 < -Ri(Pmin); $ien Eulerian relative equilibria do not exist. 

2. If (3 2 = -^i(Pmin); a unique 2-parametric family of Eulerian relative equilibria exist. 

3. If Ri(p min ) < (3 2 < 0, two 2-parametric families of Eulerian relative equilibria exist. 
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4- If (3 2 > a unique 2-parametric family of Eulerian relative equilibria exists. 

b) For k < k < and f3 2 > 0, then: 

1. If (3 2 G (£i(k), £ 2 (^))> Eulerian relative equilibria do not exist. 

2. If 13 2 = £i(k) or f3 2 = £ 2 (^)> then a unique 2-parametric family of Eulerian relative 
equilibria exists. 

3. If 13 2 > £2^) two 2-parametric families of Eulerian relative equilibria exists. 

4- If < (3 2 < £i(k) two 2-parametric families of Eulerian relative equilibria exists. 

c) For k < k < and f3 2 < 0, then: 

1. A unique 2-parametric family of Eulerian relative equilibria exist. 

d) For k < k and (3 2 > 0, then: 

1. Two 2-parametric families of Eulerian relative equilibria exists. 

e) For k < ko and (3 2 < 0, then: 

1. A unique 2-parametric family of Eulerian relative equilibria exist. 
Similarly we obtain the following result. 
Proposition 12 In the configuration S 2 S Si the following results are verified. 

a) For k > then: 

1. If (3 2 < i?i(p min ) ; then a unique 2-parametric family of Eulerian relative equilibria 
exists. 

2. If 13 2 = -Ri(Pmin); a two 2-parametric family of Eulerian relative equilibria exist. 

3. If Ri(p min ) < (3 2 < 0, three 2-parametric families of Eulerian relative equilibria 
exist. 

4- If l3 2 > a unique 2-parametric family of Eulerian relative equilibria exists. 

b) For k < and (3 2 > 0, then: 

1. If (3 2 € (£i(k), £ 2 (^)), Eulerian relative equilibria do not exist. 
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2. If (3 2 = or (3 2 = a unique 2-parametric family of Eulerian relative 
equilibria exist. 

3. If (3 2 > £2(^)5 two 2-parametric families of Eulerian relative equilibria exist. 

4- If < (3 2 < £i(k), two 2-parametric families of Eulerian relative equilibria exist. 
c) For k < and f3 2 < 0, then: 

1. If Ri(p min ) < f3 2 < 0, then two 2-parametric families of Eulerian relative equilibria 
exist. 

2. If R\(p min ) = (3 2 , then a unique 2-parametric family of Eulerian relative equilibria 
exists. 

3. If (3 2 < -^i(Pmin)) then Eulerian relative equilibria do not exist. 

4- If < (3 2 < -R2(Pmax)> then two 2-parametric families of Eulerian relative equilibria 
exist. 

5. If f3 2 = R 2 (p max ), then a unique 2-parametric family of Eulerian relative equilibria 
exists. 

6. If (3 2 > -^(Pmax)? then Eulerian relative equilibria do not exist. 

5 Stability of Eulerian relative equilibria 

The tangent flow of the equations (0) in an Eulerian relative equilibrium z e is expressed 

by 

dSz ur u 
— = il(z e )5z 

with Sz = z — z e and it(z e ) being the Jacobian matrix of (^J) in z e . 
The characteristic polynomial il(z e ) is expressed as follows 

P(A) = A 3 (A 2 + $")(A 2 + $ 2 )(A 2 + $ 2 )(A 4 + m\ 2 + n)(X s + pX 6 + gA 4 + rX 2 + s) (14) 

with $ 2 = - — j- — The coefficients present in the above polynomial are functions 

of the parameters of the problem and p, where p is taken as the root of the equation (jHJ). 
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5.1 5*2 and Si are spherical rigid bodies 

The characteristic polynomial fTlj) of il(z e ) simplifies to 

P(A) = A 5 (A 2 + $ 2 )(A 2 + $ 2 )(A 2 + $ 2 )(A 2 + K) 2 ) 2 (A 2 + p)(A 4 + gA 2 + r) 

with coefficients shown in Appendix B. 

If p > 0, q > 0, r > 0, q 2 — 4r > 0, then z e is spectrally stable. These conditions are 
not verified since r < 0. 

Proposition 13 If z e is the only relative equilibrium in the configuration SqS 2 Si of the 
zero order approximate dynamics then it is unstable. 



5.2 S2 and Si are close to a sphere 

The case in which Si (i = 1, 2) are close to a sphere will now be analyzed. In this case 
Ci — Ai w 0, and this is the reason why by applying the Implicit Function Theorem z e is 
unstable. 

If Ci — Ai is not close to zero, the coefficients of the polynomial ()14j) have very com- 
plicated expressions. Numeric calculations prove that linear stable Eulerian relative equi- 
libria exist for certain values of the parameters Ci — Ai (i — 1, 2) (see Vera and Vigueras 
[7] for details). These results are also applicable to configurations S 2 SqSi and S 2 SiS . 



6 Lagrangian relative equilibria 

6.1 Necessary conditions of existence 

If z e = (II2, HI, ITq, A e , p e x , fi e , pt) is a Lagrangian relative equilibrium, the below 
identities are verified 

A e x (V A V) e = 0, 9l I n e | 2 (A e x M e ) = (V A V) e x s 
ii e x (V M V) e = 0, g 2 I n e I 2 (A e x tf) = A e x (V M V) e 
After some calculations, it is concluded 

(A 12 ) e (A e x n e ) =0, 9l I n e | 2 (A e x n e ) = (Ai)e(A e x M e ) 
(A 21 ) e (A e x ^) =0, g 2 I ^ | 2 (A e x M e ) = (A 22 ) e (X e x M e ) 

The expressions {Aij) e are certain values of the potential function in the Lagrangian 
relative equilibrium (see Appendix A, for details). 
Concluding, the following relationships are verified 



(A 12 ) e = 0, 1 n% 

and the below proposition can be obtained. 



e ,2 (^ll)e _ (A 22 ) e 



9i 92 
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Proposition 14 Let z e = Ilf , IIq, A e , p e x , /x e , p^) be a Lagrangian relative equilib- 
rium. Then it is seen that 



92(A n ) e = gi (A22)e, (A 12 ) e = 0, \Q e \ 2 ^ {A 



ll)e 

9i 



If z e = (II2, H\, IIq, A 6 , p|, /x e , p^) is a Lagrangian relative equilibrium, then using 
the expressions of (Aj)e it is proven 

m 1 mi Pi m 2 m 2 , P2 



2 __ G(m + mi + m 2 ) 



I A e |s I /^-^A e | 3 I M e -^A e r I A e |s I M e +^A e |s I M e +^A e |* 

and 

1 n e 2 

' ' ~ Z 3 
From the previous expressions, the bellow result is deduced. 

Proposition 15 If z e = (II2, 11^, IIq, A e , p|, ^t e , p^) is a Lagrangian relative equilib- 
rium, then by denoting | A 6 |= Z, | [i e — f^|A 6 |= X, | /x e +^A e |= Y, the system of 
equations 

' X 5 - Z 3 X 2 - /^Z 3 = 



(15) 

Y 5 - Z 3 Y 2 - (3 2 Z 3 = 



has positive real solutions. 



Remark 16 Z , (3 1 and f3 2 will be the parameters that exert a strong influence on the 
study of the number of the different configurations of a Lagrangian relative equilibrium. 



6.2 Sufficient conditions of existence 

If Z has a fixed value and X, Y verify the below system of equations 

X 5 _ Z 3 X 2 _ g y3 = q 

(16) 

Y 5 - Z 3 Y 2 - (3 2 Z 3 = 

, with respect to an appropriate reference system, the coordinates of Lagrangian relative 
equilibria are completely determinated. If X = Y 7^ Z is a solution to the previous 
system, then the gyrostat S and the rigid bodies Sj (i = 1, 2) form an isosceles triangle 
when P 1 = p 2 - On the other hand, if P 1 ^ f3 2 then X ^ Y ^ Z and S t (i = 0, 1, 2) form 
a scalene triangle. 

The following proposition, whose demonstration is based on verifying the equations 
of the equilibria, shows how Lagrangian relative equilibria are when So, Si, S 2 form an 
isosceles triangle. Following the same line of reasoning, Lagrangian relative equilibria 
could be described when S , Si, S 2 form a scalene triangle. 
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Proposition 17 With respect to an appropriate reference system, z e = (H 2 , Ilf, Hq, A 6 , 
Pa> M e , p£) expressed by 

A e = (xi, yi,0), p^ = ^(-y!,xi,0), /x e = (x 2 , 2/2,0) (17) 
Pm = ^(-2/2, x 2 , 0), fig = (0, 0, wg), ng = (0, 0, C u e + 

with 

Z{m 2 -m x ) V4X 2 - Z 2 2 G(m + mi + m 2 ) 

x x =Z, 2/i = 0, x 2 = — -, 2/2 = ± ~ , K = — 

2{m\ + m 2 ) 2 z J 

are an isosceles Lagrangian relative equilibrium. The total angular momentum vector of 
the system is given by 

2 

L = (0, 0, C 2 u% + C^t + C u e + I + + yj)) 

i=i 

Proposition 18 With respect to an appropriate reference system, z e = (II|, Ilf, IIq, A 6 , 
Pl> M e > P£) 9iven by 

A e = (xi, yi,0), p^ = c?i^(-yi,xi,0), /x e = (x 2 , 2/2,0) (18) 
Pm = ^(-2/2, x 2 , 0), fig = (0, 0, u e ), U e = (0, 0, C co e + /) 

m 2 (X 2 + Z 2 - Y 2 ) - mi (Y 2 + Z 2 - X 2 ) 

xi = Z, 2/1 = 0, x 2 = — , 

2{mi + m 2 )Z 

< e ,2 G(m + mi +m 2 ) y/(Z + X + Y)(Z + X -Y){Z + Y - X){X + Y - Z) 

H)J - 7^ > 1/2 - ± ^ 

are a scalene Lagrangian relative equilibrium. The total angular momentum vector of the 
system is expressed as follows 

2 

L = (0, 0, C 2 u e 2 + C x uj\ + C u e + I + 0J e ^gi{xl + yj)) 

i=i 

Next the Lagrangian relative equilibria will be studied when S 2 and Si are spherical 
rigid bodies. 
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6.3 Lagrange relative equilibria when S2 and Si are spherical 
rigid bodies 

If 5*2 and Si are spherical rigid bodies, then C\ — A\, C 2 — A 2 . The equations (fT3j) 



are 

X 3 = Z 3 
3 _ 73 



Y A = Z 

being easily deduced from them that X = Y = Z , i.e. So, Si and S 2 form an equilateral 
triangle. It is also obtained 



n 



, 2 _ G(m + mi + m 2 ) 

" Z 3 



On the other hand, one parametrization of the former equilibria z e = (Ilf, Ilf, IIq, 
A 6 , p^, /x e , p^) is given by the relationship (|T7|) being 

Z(m 2 -mi) .VSZ 2 G(m + mi + m 2 ) 

Xi = Z, yi = 0, x 2 = — —— — , y 2 =y 2 = ±— ^— , (w ) 



2(mi +m 2 ) ' ^ a 2 ' v uy Z 3 

This parametrization of the relative equilibria will be useful to study the stability of 
the same ones. 



6.4 Lagrangian relative equilibria when S2 and Si are not spher 
ical rigid bodies 

The number of real positive solutions of the below system will be now studied 

X 5 - Z 3 X 2 - ^z 3 = 



y5 _ Z 3y2 _ = o 



(19) 



where Z, j3 1 and f3 2 are parameters. With this aim in view, let us proceed with the study 
of the number of the positive real roots of the polynomial 

p(X) = X 5 - Z 3 X 2 - (3Z 3 

according to the values of the parameters Z and /3. 

If /3 > 0, then this polynomial can only have a positive real root by applying the rule 
of the signs of Descartes. 

If (3 < 0, then we can have two positive real roots, just one positive real root or none. 
The discriminant of the polynomial, denoted by discrim(p, X), is expressed by 

12/ 



discrim(p } X) = f3Z 12 (3125(3 3 + 108Z 
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Then if discrim(p, X) < the polynomial p has two real roots, if discrim(p, x) = it 
has a positive double root, whereas if discrim(p,x) > it has no positive root. 
The discriminant is zero when the following relationship is verified 

3^20 
H 25 

From the previous results a detailed study can be carried out of the bifurcations of 
the Lagrangian relative equilibria when /3 l = P 2 . 

Proposition 19 Let z e = (II|, Ilf , IIq, A e , p^, /j, e , p^) be a Lagrangian relative equilib- 
rium and f3 l = P 2 , then 

1. If Pi > 0, a unique 2-parametric family exists, making So, Si and S 2 an isosceles 
triangle. 

2. If (3i <0 (gyrostat prolate), then: 

— 7Z 2 

al) If ^ < Pi < 0, there are two types of relative equilibria: 

• One 2-parametric family of relative equilibria, making S , Si and S 2 an 
isosceles triangle with X —Y ^ Z. 

• Two 2-parametric families of relative equilibria, making So, Si and S 2 a 
scalene triangle with X ^ Y ^ Z. 

3-^20 —1Z 2 
a2) // 25~^ 2 < ^ 1 < — 32 — ' ^ ere are ^ wo tyP es of relative equilibria: 

• Two 2-parametric families of relative equilibria, making S , Si and S 2 an 
isosceles triangle with X = Y ^ Z. 

• Four 2-parametric families of relative equilibria, making So, Si and S 2 an 
scalene triangle with X ^ Y ^ Z. 

3v // 20 

b) If Pi = ~^Z\ a unique 2-parametric family exists making S , Si and S 2 

an isosceles triangle, with X = Y ^ Z. 
3^/20 

c) If Pi < Z 2 , relative equilibria do not exist. 

On the other hand, we have that z e = (II|, Ilf, IIq, A e , p e x , /i e , p^), with equations 
(fT7j) given by 



Z{m 2 -mi) - Z 2 2 G(m + mi + m 2 ) 

xi = Z, yi = 0, x 2 = — ■ -, y 2 = ± , {u ) = — 

2 (mi + m 2 ) 2 Z 6 

describes a parametrization of the Lagrangian relative equilibria when So, Si and S 2 form 
an isosceles triangle. Similar results are obtained when So, Si and S 2 form a scalene 
triangle. 
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Proposition 20 Let z e = Ilf , IIq, A 6 , p e x , /x e , p^) be a Lagrangian relative equilib- 
rium and f3 1 7^ P 2 , then 

1. If Pi, (3 2 > a unique 3-parametric family exists, making So, S± and S2 a scalene 
triangle. 

2. If (3 1 <Q and (3 2 > then 

—7Z 2 

al) // ^ < [3 l < 0, there are two 3-parametric families of relative equilibria, 
making So, Si and S2 a scalene triangle with X 7^ Y 7^ Z. 

3^20 —7Z 2 
a2) // Z 2 < fl x < — — — , there are four 3-parametric families of relative 

equilibria, making S , Si and S 2 a scalene triangle with X 7^ Y 7^ Z. 

3v^20 

b) If fix = ^— Z 2 a unique 3-parametric family exists, making So, Si and S2 a 

scalene triangle with X = Y 7^ Z. 
3v^20 

c) If Pi < — - Z 2 , relative equilibria do not exist. 

25 

3. The case j3 2 < an d Pi > is similar to the previous one. 

4. If Pi, P 2 < then: 

—7Z 2 —7Z 2 
bl) If — — — < Pi < and — — — < P 2 < 0, there are four 3-parametric families of 

relative equilibria, making So, Si and S2 a scalene triangle with X 7^ Y 7^ Z. 

—7Z 2 3v // 20 —7Z 2 
b2) If < Pi < and — — Z 2 < P 2 < , there are six 3-parametric 

oZ Zo oZ 

families of relative equilibria, making S , Si and S 2 a scalene triangle with 
X ^Y ^ Z. 

, r, 3^20^ a -7Z 2 , 3^20" 9 n -7Z 2 , 

b3) If Z 2 < 3, < and Z l < 3 9 < , there are eiqht 

' J 25 1 32 25 2 32 ' y 

3-parametric families of relative equilibria, making S , Si and S 2 a scalene 

triangle with X 7^ Y 7^ Z. 

3^20 3^20 —7Z 2 

b4) If 3, = Z 2 and Z 2 < 9 < , there are two 3-parametric 

' J 1 25 25 2 32 ' F 

families of relative equilibria, making S , Si and S 2 a scalene triangle with 



X ^Y ^ Z. 

3/ 

m = -4 

relative equilibria do not exist. 



, c , „_ 3^20^ 3^20~ 2 . 3^20~ 2 Ja 3^20 v2 
b5) If Pi = —Z 2 andp 2 < —Z 2 or P 2 = —Z 2 and Pi < 7^Z 2 , 



19 



J. A. Vera 



Equilibria for a gyrostat in newtonian interaction with two rigid bodies 



6.5 Lagrangian relative equilibria when 5*2 and Si are close to 
spherical rigid bodies 

If Pi, (3 2 ~ 0, the solutions of (|TU|) up to the second order in the parameters f3 x and 
(3 2 are 

X = Z + ^-^- + o(Pl) 
3Z 3Z 3 KHV 



Y = Z + ^-^- + o(Pl) 



3Z 3Z 3 

Using the previous relationships, the expressions of /x e = (£2, 2/2,0) are obtained with 

(m a - m x )Z (3 X (3 2 5/3 2 5(3 2 2 2 2 

X2 ~ 2(m 2 + mi ) ~3Z + 3Z + 18Z^~ + + ^ 

V3Z .v/3 23ygg 4yg^ 23^ 2 2 

and A e = (Z, 0,0). 

7 Linear stability of the Lagrangian relative equilib- 
ria 

In the equilibrium z e the tangent flow of the equations (JTJ are 

d6z U( \x 
— = U(z e )dz 

at 

where 5z = z — z e and it(z e ) is the Jacobian matrix of (pQ) in z e . 

7.1 S2 and Si are spherical rigid bodies 

The characteristic polynomial of it(z e ) has the bellow expression 

P(A) = A 5 (A 2 + $ 2 )(A 2 + $ 2 )(A 2 + $ 2 )(A 2 + o; 2 ) 3 (A 4 + o; 2 A 2 + q) 

where 

2 G(m + m 1 +m 2 ) 27G 2 (mim + m 2 m + m\m 2 ) 

uj e = — , gr- 



and $ 2 



Z 3 ' * 4Z 6 

(C* - + I 



A 

The minimum polynomial of it(z e ) has this expression 

Q(A) = A 2 (A 2 + $ 2 )(A 2 + $ 2 )(A 2 + $ 2 )(A 2 + o; 2 )(A 4 + o; 2 A 2 + g) 
Then the bellow results are verified. 
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Proposition 21 z e is spectral stable if 

(m + m 2 + mi) 2 > 27(m 1 m + m 2 m + m 1 m 2 ) (20) 

If 

(m + m 2 + mi) 2 < 27(m 1 m + m 2 m + mim 2 ) 

then z e is unstable 

As the minimum polynomial of il(z e ) has the A = as double root, the matrix il(z e ) 
is not diagonalizable, being the following proposition verified. 

Proposition 22 The linear system 

d5z 



dt 

is unstable. 



U(z e )5z 



7.2 5*2 and Si are not spherical rigid bodies 

Similar results show that the characteristic polynomial has this expression 

P(\) = A 3 (A 2 + $ 2 )(A 2 + $ 2 )(A 2 + $ 2 )(A 2 + m)(A 2 + n)(A 8 + p\ 6 + q\ 4 + r\ 2 + s) 
The minimum polynomial of it(z e ) is expressed by 

Q(\) = A(A 2 + $ 2 )(A 2 + $ 2 )(A 2 + $ 2 )(A 2 + m)(A 2 + n)(A 8 + p\ 6 + q\ 4 + r\ 2 + s) 

The coefficients of the characteristic polynomial will be expressed in function of the 
parameters of the problem, i.e. the masses and the coefficients f3 { (i = 1,2). We conclude 
the following result by employing the Sturm Theorem. 

We conclude that the Lagrangian relative equilibria are spectral stable (lineally stable) 
if the following conditions are verified 

p 2 q 2 - 3rp 3 - Qp 2 s - 4q 3 + Upqr + lQqs - 18r 2 > (> 0) 
p 2 qr — 48sr — 9sp 3 + 32pqs — Aq 2 r + 3pr 2 > (> 0) 
r, s > (> 0), 3p 2 -8q>0(> 0), pr - 16s > (> 0) 
m, n > (> 0), discrim(h) > (> 0) 

If Si (i — 1,2) are arbitrary rigid bodies, the former conditions are very complicated 
expressions in the parameters of the problem and can only be studied by means of nu- 
merical analysis. 
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If Si (i = 1, 2) are close to a sphere, the coefficients of P, up to the first order in the 
parameters (3 X and (3 2 are expressed as follows 

G(m + m 1 + m 2 ) G{m + m l + m 2 ) 
m= ^ +o(/3 1 ) + o(J3 2 ), n = — + o{p x ) + o(/3 2 ) 

27G 3 (m + mi + m 2 )(mim + m 2 m + mim 2 ) , . , . 
r = ^ + + o(/3 2 ) 

G 2 (4m„ + 4m? + 4m? + 35m m! + 35m m 2 + 35mim 2 ) . . . . 
9 = + o(/3i) + o(P 2 ) 



If the function 

81G 4 m (m + mi + m 2 ) 2 (/3 1 mi + /3 2 m) 
S= 4^ 

is positive and 



+ o(p 2 1+ Pl) 



(mo + m 2 + mi) 2 > Tl{m\m§ + m 2 mo + mim 2 ) 
z e is linearly stable. Then if 5j (z = 1, 2) is close to a sphere and 

m x {Ci - A x ) + m 2 (C 2 - A 2 ) > 

z e is linearly stable when pOj) is verified. If mi(Ci — Ai) +m 2 (C 2 — A 2 ) = 0, second order 
terms of the coefficients in (3 1 and f3 2 are necessary for the study of the problem. 



8 Conclusions and future works 

In this paper it has been investigated some important periodic solutions of the dy- 
namics of a gyrostat in Newtonian interaction with two symmetric rigid bodies. With the 
hypotheses formulated in the introduction of this paper, working in the double reduced 
space of configuration of the problem, both the equations of motion and the those which 
determine the relative equilibria have been derived. 

Two families of relative equilibria, Eulerian and Lagrangian are studied. The Eulerian 
relative equilibria have been completely determinated by a polynomial equation of degree 
nine. The bifurcations of these equilibria have been carried out when mo is very small. 
The bifurcations of the Lagrangian relative equilibria have been completely investigated 
and expressions for the Lagrangian relative equilibria, when both solids are close to a 
sphere, are provided. These expressions are useful for the posterior study of the stability 
of the same ones. 

The instability of Eulerian relative has been proven and necessary and sufficient con- 
ditions for lineal stability of Lagrangian relative equilibria are provided. 
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Different results, which had been obtained by means of classic methods in previous 
works, have been generalized in a different way, and other results, not previously consid- 
ered, have been studied. Some interesting numeric results have been detailed in Appendix 
B. 

Numerous problems are open, and among them the study of the "inclined" relative 
equilibria, in which form an angle a^O, n/2 with the vector A e x /x e , is considered. 

The methods used in this work are susceptible of being used in similar problems. The 
nonlinear stability of the relative equilibria obtained here is the logical continuation of 
this work. 



9 The (A,-)e, i,j = l,2. 

The bellow expressions of the potential V are obtained. 

Gmim 2 \ e Gm m 2 ^. a l (^-J^^) Gm mi^ a 2 (/i e +f^A 6 ) 



^ V \ V ) e — I \e Is i\/r / j\ ..p mo\e ' 



A e I 3 M 2 ^1 /x e -^A e p+ 3 M 2 ^ I /x e +^A e | 2i + 3 

I i=0 -^2 ' Z i=0 1 ^ M 2 1 

(V M V) e = Gm 2^1 ,,e_I22 A e |2i+3 + 2^1 ,.e T \e | 2i+3 
\i=0 1 M I i= 1 ^ + M 2 A I 

where = a 2 = m , a\= f3 l = 3/2(Ci - A x ) and a 2 = /3 2 = 3/2(C 2 - A 2 ). 
The following identities are also verified 

(V A V) e = (A n ) e A e +(A 12 ) e /x e , (V M V) e = (A 2l ) e X e +(A 22 ) ef i e 

where 

Gm 1 m 2 Gm m 2 f a i \ Gm ml ( a 2 

( ll)e = ~wf + ~w \h'\ » e -M xe \ 2i+ v + ~ JI f~ \h'\ » e+ ^ xe i 2i+3 . 

_ Gitlqitli / a? \ Gm m 2 ( ^ a] 

' 12 ^ e ~~ A/L 2^1 „e i 22L\ e |2i+3 JLf„ 2^1 „e_: 



M 2 1^1 /x e +^A e | 2l + 3 / M 2 1^1 /x e -^A e | 2 *+ 3 

^ \i=0 1 ^ M 'i 1 / ^ \i=0 1 ^ M 2 

(A 12 ) e = (A 21 ) t 



1 



/ ■ a 1 a 2 

(A 22 ) e = Gm J^l n e_rnl X e , 2i+3 + I ,.e 7 mAe , 2i+3 



i=0 1 ^ M 2 I j=o I " M 2 ' 
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10 Some numerical results 

In order to obtain the values of Cj — A-i, (i = 1,2), the below relationships will be 
utilized 

C 1 -A 1 = (l-v) (^)V 

where es i and ps { represent the equatorial and polar radius of Si (i = 1,2), jf' = |£j, 

respectively, and Si = — 1 -. Si, S2 are considered to be homogeneous ellipsoids. The 

distances are measured in kilometers. 



^^^(mo -> 0) 


Without oblat. 


Oblat. of S 2 


Oblat. of S 2 


and Si 


Earth-Moon-S'o 


448879.206 


448879.221 


448879.251 




Mars-Phobos-^o 


9414.945 


9414.958 


9414.958 




SoS 2 Si(m -> 0) 


Without oblat. 


Oblat. of S 2 


Oblat. of S 2 


and Si 


^o-Earth-Moon 


381679.691 


381679.763 


381679.763 




S'o-Mars-Phobos 


9310.642 


9310.666 


9310.668 




S 2 S Si(mo -> 0) 


Without oblat. 


Oblat. of S 2 


Oblat. of S 2 


and Si 


Earth-So- Moon 


326409.744 


326409.780 


326409.751 




Mars-S -Phobos 


9339.156 


9339.196 


9339.196 





24 



J. A. Vera 



Equilibria for a gyrostat in newtonian interaction with two rigid bodies 



11 Coefficients of the characteristic polynomial in 
Eulerian relative equilibria S0S2S1 

The coefficients of the characteristic polynomial (fT4*j) are 

2 _ G ((m 2 + mi) p 4 + (2 mi + 2 m 2 ) p 3 + (m 2 + mi) p 2 - 2 m p - m ) 

_ G((m 2 + 4m + mi)p 3 + (3m 2 + Qm )p 2 + (4m + 3m 2 )p + m + m 2 ) 
P " (l + p) 3 P 3 A 3 

G((—2mip 4 m 2 + (— Im^mi + + m| — 2mim 2 — 2m§m 2 )p 3 

q= (OTpJM) + 

(3m 2 + mim 2 — 6momi)p 2 + (—mim 2 + 3m 2 + 2m m 2 — 4m m 1 )p 

((1 + P) 3 P 3 A 3 ) + 
m 2 — m^mi + mom 2 — mim 2 )) 
((l + p) 3 P 3 A 3 ) 

_ G 2 (aip 4 + a 2 p 4 + ct3P 2 + a 4P + a 5 ) 
11.1 Coefficients a$ (i = 1, . . . , 5) 



T R O R R O 

ai = — 42m 2 mi — 48m 2 mo — 147m 2 m 1 — 336m 2 mimo — 129m 2 m 
— 207m 2 m 3 — 782m 2 m^m — 673m 2 mimo — 81m 2 mQ — 150m 2 m 4 
— 869m 2 m^m — \325m\m\m 2 ] — 378m 2 m 1 mo — 64m 2 m^ 
— SlSmgm^^mo — 1270m 2 m 1 m — 702m 2 m 1 m — 14m 2 m 1 
— 165m 2 m^m — 610m 2 m 4 mQ — 648m 2 m^mQ — 24m 2 m^m 
—119m 2 mlml — 297m 2 m 4 m 3 ) + 2m t lm 2 ] — 54m ! j'm 3 ) 
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a 2 = —QOm^mx — 54m 2 m — 2A3m 2 m 1 — 474m 2 m 1 m — 173m 2 m 
-399m 5 2 ml - 1345m 2 mim - 999m^mimg - 135m 5 2 ml - 32%m\m\ 
— \SAQm\m\mQ — 2223m 2 m\m^ — 648m 2 m 1 m ! — 138m 2 mi 
-l3QAmlm\m - 2hQQm\m\ml - \2A2rr? 2 m{m\ - 2Am 2 2 m\ 

OK O/IO O Q Q f\ 

— 536m 2 m 1 m — 1530m 2 m 1 m — 1188m 2 m 1 m — 90m 2 m 1 m 
— 477 m^m^m^ — 567m 2 m^mQ — 56m ( j i m ! — 108m^mQ 

03 = — 42m 2 mi — 36m 2 m — 183m 2 m 1 — 342m 2 mimo — 93m 2 m 
—3A9m\m\ — 1097m 2 m^m — 630m 2 m 1 m ! — 81m 2 mQ — 358m 2 m^ 

A Q /IOO /I Q Q ^ 

— 1776m 2 m 1 mo — 166m 2 m 1 m — 405m 2 mim — 189m 2 m 1 
— 1614m 2 m 1 mo — 2256m 2 m 1 m — 810m 2 m 1 m — 31m 2 m 1 

Qpr 0/10 OQQ 7 

— 827m 2 m°m — 1683m 2 m 1 m — 810m 2 m^mQ — 6m 2 m( 
— 228m 2 m 1 m — 666m 2 m 1 m — 405m 2 m 1 m — 30m 1 mo 
-81m^m^ - lllmfmjj 

a 4 = — 12m 2 mi — 12m 2 m — 5Qm 2 m 1 — 114m 2 mim — 24m 2 m 
— 130m 2 m^ — 387m 2 m ; j ! mo — 162m 2 mimQ — 179m 2 m\ 
— 687m 2 m^m — A32m\m\m 1 {i — 140m 2 m^ — 693m 2 m^m 

00Q OR O ^ 0/10 V 

— 588m 2 m 1 m — 52772277^ — 387m 2 m 1 m — 432m 2 m 1 m — 6m 2 m 1 
— 108m 2 m^m — \Q2m 2 m\m^ — 12m]m — 2Am\m^ 

05 = — (mo + m 2 )(18mom 2 + 12mim 2 + 94m 2 momi + 36m^m 2 

/10 /lOi/l 4 ? Q Q 

+81m 2 m m 1 + 168m 2 m m 1 + A2m 2 m\ + 128m 2 m m i 
+27m 2 m 1 + 15m 2 mf + 31m 2 m m^ + 12Qm\m^m\ + 18mQm 2 
+54m 2 mgm^ + 12m 2 m m^ + hm 2 m\ + 7m^m + 9m ! m ! j i 
+144m 2 m ! m : j ! ) 
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